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Abstract 

An effective model is suggested for an asymmetric spin ladder with dimerized 
rungs. Magnon mode truncation originated from magnon decay (recently observed 
in the ID compound IPA — CUCI3) is naturally described within this model. Using 
Bethe Ansatze we studied a one-magnon sector and obtained relations between 
interaction constants of the model and experimentally observable quantities such as 
the gap and truncation energies, spin velocity and truncation wave vector. It is also 
shown that the structure factor turns to zero at the truncation point. 

1 Introduction 

A spin ladder with strong antiferromagnetic rung coupling gives an ideal example of 
a gapped spin-dimerized system pQ. Really the majority of spins in its ground state 
are coupled in rung-singlets (rung-dimers) so that the relative coupling energy estimates 
a value of the gap. By this reason all low-temperature effects depend on dynamical 
properties of states with a few number of excited rungs. Theoretical study of excitations 
in spin ladders with strong antiferromagnetic rung coupling was developed in a number 
of papers [2]- [5]. It was pointed out that the lowest excitations form a coherent magnon 
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branch. When the gap energy is smaller than the energy width of the magnon zone the 
latter may intersect the two-magnon scattering continuum. For a symmetric spin ladder 
(with equal couplings along both legs as well as along both diagonals) these two sectors 
do not hybridize so a one-magnon state is always stable. The situation is quite different 
for an asymmetric spin ladder with non equal couplings along legs or along diagonals. As 
it was pointed in [6]- [9] the coupling asymmetry entails hybridization between the "bare" 
(related to a symmetric case) one- and two-magnon sectors. If the system has a wide 
magnon band intersecting with the two-magnon scattering continuum this hybridization 
results to magnon instability and truncation of the magnon mode at some value k trunc of 
wave vector. Experimentally a magnon mode truncation was observed in neuron scattering 
from ID compound IPA - CuCl 3 ((CH 3 ) 2 CHNH 3 CuCl 3 ) [6J. The latter is considered as 
an asymmetric spin ladder with strong ferromagnetic rungs and is effectively equivalent 
to a ID Haldane antiferromagnet. 

In this paper starting from an asymmetric rung-dimerized spin ladder we present an 
effective model which produce an explicit realization of magnon mode truncation related to 
magnon decay. Within our model we study the one-magnon excitations and obtain explicit 
relations between the coupling constants and experimentally observable quantities. 



2 Hamiltonian for an asymmetric spin ladder 

In the present paper we shall study an asymmetric spin ladder with the following Hamil- 
tonian 

H = y^#n,n+l) (1) 

n 

where # n , n+1 = + H l * 9 n+1 + H*™ 9 +1 + H° v ° +1 + #™ and 

J_i 
2 



Tjrung _ -L o _j_ O Q ^ 

n n,n+\ ~ "7TV°l,n ' °2,n T Oi jn+ i • 02,n+l, 



Hnfn+l = J\\{^>l,n ■ Si jn+ i + S 2 , ra ' S 2 ,n+l), 
^n"n+\ = JdSl,n ' S 2 ,n+1, 

^n,n+l = ^c((Sl, n ■ Si jn+ i)(S 2jn • S 2;rl+ i ) + (s l,n ' S 2>n ) 

x (Si, n +i ■ 82,71+1) — (Si >n • S 2jn+ i)(S 2in • Si )Tl +i)). (2) 

Here Sj tU (j = 1,2) are the S — 1/2 spin operators related to n-th rung. The auxiliary 
term H norm = J norm I (I is an identity matrix) is added for the zero normalization of the 
ground state energy. 
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The following condition, 

J d +J c = 2J h (3) 

suggested in [4], guarantees that the vector, |0) n (g> |0) n+ i (where |0) n is the n-th rung- 
singlet or equivalently rung-dimer) is an eigenstate for H n<n+ i, so the vector 

io>=ni°>»> w 

n 

is an eigenstate for H. An additional system of inequalities, 

5 3 

J± > 2J||, J_L>-J C , J± + J\\>-^Jc, 



3J ± -2J||-J C > J.J\ -4Jj.J|| +20J,f - 16J||J C + 4J C 2 , (5) 



together with a condition J norm = 3/4Jj_ — 9/16J c , guarantee that the vector (4) is the 
(zero energy) ground state for H. The full system of the "ground state tuning" conditions 
(3), (5) belongs to the mathematical basis of our model. 

For J d = the Hamiltonian H commutes with the operator, Q = \ ^ n (Si in + S2, n ) 2 , 
considered as a number of "bare" magnons [5], [9]. Therefore the Hilbert space splits on 
an infinite sum: Ti. = J2m=o^- m ^ wnere Q\n m = m - The subspace H° is one-dimensional 
and generated by the ground state (4). 



3 Spectral problem for the reduced Hamiltonian re- 
lated to the effective model 

Despite the ground state (4) for the Hamiltonian (l)-(3), (5) is known it is not clear how 
to obtain its excitations. In the symmetric case [4j,[5] the one-magnon state corresponds 
to H 1 but even for a small asymmetry it already lies in *Y^ = $'H 2n+1 |9J. By this reason 
the related spectral problem seems to be unsolvable. However for a strong rung coupling 
the states with rather big number of "bare" magnons have a large energy and therefore 
may be effectively reduced. In the first order with respect to the dimerization energy the 
reduced Hilbert space 

H red = n o QH i (sn 2 contains additionally to the ground state (4), 
only the "bare" one- and two-magnon sectors. The corresponding effective Hamiltonian 
H e ff is defined as the restriction of H on Ji red or, 

fteff = p(0,l,2)£p(0,l,2) ; (6) 
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where p(°' 1,2 ) is the projector on H red . 

A general S = 1 excited state for H e ^ related to a wave vector k and the energy E(k) 
is superposition of "bare" one- and two-magnon components, 

|*)« = -^X;[a(A^^ (7) 

where = (S° — S"^) |0) ra and means an infinite product of dimers related to the 
remaining rungs. The normalization factor Z(k) is defined as, 

oo 

Z\k) = \a(k)\ 2 + 2Y,\b(k,n)\ 2 . (8) 

The system of Shrodinger equations on the amplitudes a(k) and b(k, n) directly follows 
from the local action of the operator H n<n+ i, 

H n , n+1 \o) n \ir n+1 = (ij ± - ^J c )|0) n |l)^ +1 + ^|l)^|0) n+1 - ^e a M\l)Z +1 , 

H n:n+1 \l)%\0) n+1 = (Ij ± -^J c )|l)«|0) n+1 + ^|0) n |l)« +1 , 
H n , n+1 e aPl \ 1)3 = (J ± - J,,- J c /4)£^ 7 |l)^|l)^ 1 + iJ d |0) n |l)« +1 . (9) 

From (7) and (9) one can obtain an infinite set of recurrent equations, 

k 

(2J ± -3J c )b(k,n) + J c cos-[b(k,n-l) + b(k,n + l)]=E(k)b(k,n), n > 1, (10) 

related to non- neighbor excited rungs and two additional equations related to neighbor 
rungs, 

3 k 
(J± J c + J c cos k)a(k) + iJd cos — b(k, 1) = E(k)a(k), 

(2J J .-^J c -^)6(fc,l) + J c cos^(fc,2)-^oos^a(fc) = E(fc)6(fc, 1). (11) 

For a coherent excitation originated from the hybridization of the one-magnon and 
bound two-magnon states there must be, 

lim b(k, n) = 0. (12) 

With regard to this condition the Eq. (10) has the following general solution, 

b{k,n) = B{k)z n {k), (13) 
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where, 



and 



\z(k)\ < 1. 



E{k) = 2J ± - 3J C + J c (z{k) + -L) 
From (14) and (15) follows that, 



(k)J C °4 



(14) 
(15) 



lmz(£;) = 0. (16) 
Substituting (13) and (15) into (11) we obtain a pair of equations on a(k) and B{k) 



represented in the following matrix form, 



M(k) 



a(k) 
B{k) 



0, 



(17) 



where 



M(k) 



( 3 

2 

V 



Jc 



J c cos k — J c (z(k) + ^pjy j c °s ^ — J. 



k 



iz{k)Jd cos 



k 



(jJc-if)z(k)- J c COs| 



(18) 



The Eq. (17) is solvable only for detM(k) = 0, or, 



k 

z 2 (k)J c cos — + ( J 



2 



x 



J c cos A; 



(| J c- Jd)z(k) -2 J c cos | + z 2 {k)Jl cos 2 | = 



)*(*) 



+ J c cos — 
2J 



(19) 



The Eq. (19) added by the conditions (14) and (16) completely defines the coherent 
spectrum for H e ^ . The truncation originates from a failure of any of the conditions (14) 
or (16). For the first possibility the truncation wave vector k trunc coincides with the, 
critical wave vector k c defined as, 

|*(fc c )| = l. (2°) 

For the second one it coincides with the branching wave vector kb related to passing of 
solutions of the Eq. (19) into the complex plane. 

In order to clear the nature of the truncation point for an arbitrary set of coupling 
parameters (however limited by (3) and (5)) let us first examine the case when the con- 
dition (20) is satisfied just at the branching point. In other words we are interesting in 
k c = h = he when the Eq. (19) has a twice-degenerate solution z{kb c ) so that the same 
one has the equation obtained from (19) by differentiating of its left side with respect to 
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z(k). Using an auxiliary variable / = z(kb c ) cos and taking in mind that according to 
(16) and (20) z 2 {kb c ) = 1, we represent (at k = kb c ) the Eq. (19) and its "derivative" 
equation as the following system, 

+ (J 2 d + Wc - 7J 2 )f 2 + J C (AJ C - 2J d - 2J ± )f + (p c 
2J 2 f + (J 2 + 2J d J c - 5J 2 )f 2 + J c ( 7 -J c - 2J d - J ± )f + (p c 
which is solvable only for, 

J c Jd(2J±- J c )(3J c -2J d ) = 0. 

(The left side of (22) was obtained from the resultant of the two polynomials in the left 
sides of (21)). The solution J c = of the Eq. (22) is not interesting because in this case 
the Eq. (19) is singular and solvable only for k = n. The solution J c = 2Jj_ is inconsistent 
with (5). The solution 3 J c = 2J d is artificial because in this case cos | factorizes from the 
left side of (19), and therefore at k = % the Eq. (19) is identically satisfied for all z(tt). 
The solution J d = relates to zero asymmetry when the corresponding truncation wave 
vector k , 

-4-Wx- 1 ). (23) 

may be easily obtained from the Eq. (19) which also gives, 

z(k ) = -1. (24) 

The formula (23) has a clear physical interpretation. Really as it follows from the results 
of the Refs. [1] and [5] (related to symmetric spin ladders) at k = k$ the "bare" one- 
magnon branch with dispersion E™^ 9 e n {k) = J± — |J C + J c cosk intersects the lower bound 
of the scattering two magnon continuum [5], 

E ^a 9 n,lo W{k) = 2 J± _ 3 ^ _ ^ k_ (25) 

The above result confirm the general statement suggested in [0]-[E] that even an extremely 
small asymmetry may change drastically a magnon mode. As it follows from (23) at 
J d — » the truncation occurs only for 9J C > 2J^. Since J± > the parameter J c also 
must be positive. 

In order to find a nature of the truncation at J d ^ let us study an evolution of z(kb) 
for small J d . If the condition (14) is satisfied for k = kb then the truncation originates 



-Jd)(J±- 



Jr 



0, (21) 



(22) 
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from branching and k trunc = k b . However in the opposite side for |z(A;&)| > 1 it will be 

ktrunc k c . 

Taking for Jd/ J c 1 and k ko the following infinitesimal representation z(k) = 
— 1 + e(k), using the following notations t(k) = cos/c/2, t = cosA;o/2 and the formula, 

J 3 

-f- --cos k = 2(t 2 + t -t 2 (k)), (26) 

which follows from (23) we obtain from (19) by omitting the term e 3 (k) the following 
equation, 

a{k)e 2 {k)+(3{k)e(k)+ 1 {k)=0. (27) 

Here 

_ 1 , t(fc) o (t(fc)-t )(t(fc)+t + l) Jjt(fc) 
„ m _ o (t -t(fc))(t(/c)+t + l) /o J, 2 t(A:) 

w " 2 ^) (2+ a7 )+ jfA7' 

7(A;) = 2 (1 + a7 } " 2j|aT' (28) 

and Ai = 3/4- J d /(2J C ). 

The branching wave vector is characterized by the following condition, 

D(k b ) = (3\k b ) - 4a(hh(h) = 0. (29) 

After its linearization with respect to small parameters t(k) — 1 and Jj/ J 2 this equation 
reduces at first to, j(kb) = and then to, 

-i^+ imM^r (30) 

According to (27) and (29), e(k b ) = -/3(k b ) / '(2a(k b )), or using (30) and (28), 

Since e(k b ) < the condition (14) fails for z(k b ). Therefore for Jj <C J 2 there must 

be, 

ktrunc = k c . (32) 

Since for J d ^ the wave vector k c evolve continuously from k the Eqs. (24), (20) and 
(16) give, 

z(k c ) = -I. (33) 



7 



Despite the Eqs. (32) and (33) were proved for Jj <C J% they are right for all Jd- Really 
if for some region of Jd it will be k trU nc = kb then there must be a point where k c = kj,. 
But as it was shown above ko is the only one point of such type. 

The Eqs. (15) and (33) give the following representation for the magnon energy at 
the truncation point, 

k 

Efrunc = 2 J± - 3 J c - 2J C COS y . (34) 

The magnon branch approaches the bottom of the two-magnon continuum tangen- 
tially, 



d_ 

dk' 







k=k c dk 



E 2raa g n,lo W{k) = " E (h) =J c shl^. (35) 



The Eq. (35) may be easily derived from (25) and (15) using an auxiliary relation, 

+ = ' ( 36 ) 

ok \ z(k)J k=k c 

which follows from (33). The same result was obtained in [7J by a different approach. 

Let us notice that the singularity at to+A\ = in the formulas (30) and (31) originates 
from a resonance between the one-magnon and bound two-magnon states. Really for 
t = — Ai the Eq. (19) has the thrice-degenerated solution related to both these states. 
This special case is not considered in the present paper. 



4 Magnon dispersion near the gap 

Let us turn to the opposite side of the spectrum related to k = ir. As it follows from (19) 
z(k) is an odd function and 2(71*) = 0. Therefore for k w tt we may put, 

z{k) ^ Zi {tt - k) + z 3 {n ~ k) 3 . (37) 

Then from (15) and (37) follows that for z\ — z±/12 — Z3 > the dispersion at k ~ it takes 
the form, 

E(k)*E gap (l + ^(n-ky), (38) 
where the gap energy E gap and the spin velocity v sp i n are given by 



rp _ 9 t _oT_i_^ c v spin 

' gap 



2zi E, 



1 _ 

1 



\ 2J ± -3J C + ^ 
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Since both E gap and v sp i n may be obtained by an experiment [6] we shall express them 
explicitly from the coupling constants. 

Substituting (37) into (19) we obtain the following system of equations on the coeffi- 
cients Z\ and z 3 , 



J i 



Jc 

' 2\ 



Jr 



J 



+ 



Jr. 



c / _2 



1 



z i- — ) +z 3 (J±- — 



2 V 1 12 
The Eq. (40) has two solutions, 

magn _ 
Z l — 



Jc 



~J C — Jd }Zi — J c 



;J C ~ Jd)zi - J< 

Jc 

12J 

T 2 z 2 

4 



0. 



(40) 



2^c - Jd )z 3 + — 



+ 



Jr 



2J, - J 



bound 



2J r 



0. (41) 



(42) 



3 J c — 2J^ 

related to magnon and bound two-magnon branches [I], [3]. According to the first equation 
in (42) zf aan {J ± - JJ2) = ~J c z 3 nagn /(2z 7 [ iagn ) so from (41) follows, 



magn 



y magn 



72 



and according to (39) one can obtain, 



J C (4J± + J C - 2Jd) 



(43) 



E, 



gap 



J I — — Jr 



J spm 



2J r 



E, 



<i«i> 



2 J, - 5 J c 



J 2 



J C (4J ± + J c - 2J d 



(44) 



As it follows from (44) the point k = n corresponds to an energy minimum (the gap) only 
for J c (4Jj_ + J c — 2J d ) > J 2 . (According to the comment after the Eq. (24) we suppose 
that J c > 0.) 

Using (33) and (34) we may represent the Eq. (19) in the point k = k c as follows, 



(3 



trunc 



k 

2E gap ) cos 2 y 



1-cos- 



+ E 



trunc 



2E, 



gap, 



3 + 4 cos • 



l( 1 - cos — 

2 



gap 



~Jd 



E 



trunc 



Jjcos 2 f. 



(45) 



where the parameters J± and J c are excluded by (34) and (44). The Eq. (45) may be 
used for obtaining the parameter J d directly from an experimental data. 
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5 One-magnon dynamical structure factor near the 



threshold 

We use the following representation for the dynamical structure factor (DSF), 

Sapfau) = hm ^(OlS^^I^^IS^-q)^)^-^). (46) 



Here S(q) is the spin Fourier transformation associated with the two dimensional vector 
q = (q, q rU ng) with leg and rung components. Since the latter has only two possible values 
and 7r we may study them separately, 

S( 9 ,0) = ]Te-^(S lin + S 2 , n ), S(q,n) = J^V^S^ - S 2 , n ). (47) 

n n 

According to the following pair of relations, S(g,0)|0) = 0, S(— q,ir)\0) = X] n e ' 9 " '■■■\^)n---, 
we may reduce the matrix elements in (46) 



(/i|S(g,0)|0) = 0, «(k\SP{-q^)\0) = S aP S kq ^^-, (48) 

therefore, the DSF has purely diagonal form, S a p(q,7r,u>) = 8 a pS(q,ir,uj), while the one- 
magnon contribution is purely coherent, 

S magn (q,n,u) = A magn (q)5(u - E ma ° n (q)), (49) 

where 



A-magn ((?) 

According to (8) and (13), 



"rnagn 



Zmagn ((?) 



(50) 



v majn(^) \j \^magni}^)\ Z^(k^) ' ) 

For q — > k c , it will be A magn (q) oc 1 — z 2 (q), so as it follows from (33) A magn (k c ) = 0. The 
same result was obtained in [7j by different approach. 

Finitely let us notice that a rather similar effect of hybridization between magnon and 
phonon modes was studied in [10]. However in the latter case a magnon mode does not 
truncate (because there is no decay) and therefore the corresponding structure factor does 
not turn to zero. 
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6 Summary and discussion 



In this paper for a rung-dimerized asymmetric spin ladder we suggested an effective model 
which neglects all states with n > 2 bare magnons. Using Bethe Ansatze we studied the 
effect of magnon mode truncation resulting from magnon decay and clarified its mathe- 
matical nature (see the Eq. (32)). We obtained the four equations (see (34), (44) and 
(45)) coupling the interactions constants of our model (namely J±, J c and J^) with the 
truncation wave vector, gap and truncation energies and spin velocity. 

Of course the neglect of the states with n > 2 bare magnons is a rather rough approx- 
imation. Really an intersection between the one- and two-magnon scattering modes is 
possible only for a wide band system. In this case the bare n > 2 zones also lie not so far 
from the magnon mode and therefore give a rather essential contribution to it. However 
if we concern only on the gap and truncation points then our model produces a good 
approximation. Really near the gap the magnon energy is minimal and lies far below the 
bare n > 2 magnon modes. For example as it follows from (44) for Jd -C J±, J c even the 
n = 2 correction is small. From the other side since the 1 — > 2 decay threshold lies on a 
finite distance below the 1 — > 3 one the latter is not sufficient at the vicinity of the trunca- 
tion point where the parameters (E(k c ) — E(k))J E gap and (E(k c ) — E{k))/ (E trU nc — E gap ) 
are small. Therefore the infinitesimal analysis of the Sect. 3 (Eqs. (27)-(32)) gives the 
right picture of the truncation (the Eq. (33)). 

There is only one known asymmetric rung-dimerized spin ladder compound namely 
the CuHpCl (see [H] and references therein). However the effect of truncation was not 
observed in this material. This fact is clear because the gap energy in CuHpCl (0.9 mev) 
is bigger than the magnon bandwidth (0.5 mev) so the magnon mode does not intersect 
with the scattering two magnon continuum. 

Despite none wide-band asymmetric rung-dimerized spin ladder compound was found 
up to now we suppose that this may likely happen in not so remote future. Then the 
results of our paper probably will be useful for a theoretical study of such compound. 

The author is grateful to S. L. Ginzburg, S. V. Maleyev and A. V. Syromyatnikov for 
the interest and helpful discussion. 
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